Influence of confinement on granular penetration by impact 
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We study experimentally the influence of confinement on the penetration depth of impacting 
spheres into a granular medium contained in a finite cylindrical vessel. The presence of close lateral 
walls reduces the penetration depth, and the characteristic distance for these lateral wall effects is 
found to be of the order of one sphere diameter. The infiuence of the bottom wall is found to have 
a much shorter range. 
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The motion of a solid sphere falling under its own 
weight into an unbounded viscous fluid has been known 
for more than one century to be characterized by the lim- 
iting Stokes velocity, and the influence of a confinement 
by close walls on the sphere velocity is also known for 
many years both theoretically Ikl^, 0] (see Ref. Q for 
a review) and experimentally [5|, l6|. This wall influence, 
due to the long range of hydrodynamic forces for low 
Reynolds number flows, results in a larger fluid friction 
force on the sphere. The velocity of a sedimenting sphere 
is thus known to be reduced near a wall either parallel or 
perpendicular to the sphere motion by a correction factor 
that is linear with d/l at low d/l {d/l < 1), where d is 
the sphere diameter and I is the distance from the wall. 
In all these cases, the well known equations of classical 
fluid mechanics serve as a reliable support to understand 
and predict the corresponding motions with a fluid fric- 
tion always proportional to velocity at small Reynolds 
numbers. 

Considering now the motion of a falling sphere into 
granular matter, the existence of a non zero solid fric- 
tion force at vanishing velocity explains that the sphere 
driven by its own weight cannot reach a constant limit ve- 
locity but stops suddenly at a given depth. As the precise 
rheology of granular matter is far from being well under- 
stood and does not benefit from well-accepted equations 
despite recent important progress 0, the predictions for 
the stopping distance are far from easy. Many studies 
have been performed on the penetration depth 5 of an 
impacting sphere into a granular medium. It emerges 
from those studies [1, 0, that i5 follows a power law of 
the form 5/d oc [p/ PgY {H / d)"^ , where Pg is the density 
of the grains, d and p are respectively the diameter and 
the density of the sphere and H the total falling distance 
covered from release to rest. The first power exponent is 
/3 ~ 0.5 whereas the second, a, varies between 0.3 and 
0.5 depending mainly on the range of impact velocities 
and slightly on the packing fraction of the layer 0, [^, • 
Indeed, as the packing varies from dense to loose random 
packing, the packing density does not vary so much but 
the contact network and force network do vary signifi- 
cantly. 

The knowledge and the expression of the forces exerted 



by the granular medium, and responsible for the sphere 
deceleration and stop, are still the subject of an intense 
debate [ll|, [l^] , with different terms of frictional and col- 
lisional origins. Recently, the effect of a solid wall normal 
to the sphere motion was shown to account for an expo- 
nential increase of the force on the sphere in the close 
vicinity of the wall p^[ll|. 

To our knowledge, no study deals in detail with the in- 
fluence of close walls on the granular penetration depth, 
except first results of Ref. [15] . In the present paper, after 
checking the usual "unbounded" case, we investigate ex- 
perimentally the influence of the walls, either parallel or 
perpendicular to the sphere motion, for a sphere dropped 
onto a granular layer. In some situations, a coupling be- 
tween the effects of parallel and normal walls may exist, 
but to understand their respective roles, each is studied 
separately in this paper. 



I. EXPERIMENTAL SETUP 

The experimental setup is sketched in Fig. [TJ The pen- 
etration depth 5 is investigated by dropping a spherical 
projectile of diameter d and density p onto a fine granu- 
lar material in a cylindrical container of diameter D and 
height h. The granular medium consists of glass beads 
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FIG. 1: Sketch of the experimental setup and notations. 
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(density pg ~ 2.5 x 10^ kg m~^) slightly polydisperse in 
size, with a diameter range of 300-400 /im. Before each 
drop, the granular medium is prepared by gently stir- 
ring the grains with a thin rod. The container is then 
overfilled and the surface levelled using a straightedge. 
We have checked that this preparation leads to repro- 
ducible results with only small variations. The grain size 
is much smaller than the falling sphere diameter d so that 
the granular medium can be considered as a continuum 
medium. 

Different sphere materials and sizes have been used to 
bring out the influence of the sphere density p and of 
the sphere diameter d on the penetration depth 5. Steel 
spheres are initially maintained by a magnet at a distance 
h above the granular surface. Non-metallic spheres are 
held by creating locally a vacuum at the top of the sphere. 
Both apparatus allow one to drop the spheres without 
any initial velocity nor spinning motion. The sphere is 
released directly above the center of the container and 
fall along the container axis. The penetration depth 5 
is then measured with a thin probe by locating the top 
of the sphere with a precision better than 1 mm. The 
impact speed is tuned by varying the drop height h from 
to 0.5 m; the corresponding velocity v at impact, given 
by t; = "v/SgTi where g is the gravitational acceleration, 
thus ranges from to 3 m s~^. In the following, the 
results will be analyzed using the total falling distance 
H ^h + 6 (Fig.[T]). 
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FIG. 2; Normalized penetration depth 5/d as a function of 
the normalized total falling distance H /d. Experimental data 
for steel spheres (p ~ 7.8 x 10'^ kg m~'^) of different diameters: 
(V) d = 5 mm, (o) d = 10 mm, (□) d = 19 mm, (o) d = 20 mm 
and (A) d = 40 mm. (— ) Power law fit {5/d) oc {H/d)" 
with a — 0.4. The shaded region corresponding to H < 
S is by definition not allowed. Inset: {5 / d){H / d)~^''^ as a 
function of the density ratio p/ pg for spheres of diameter d = 
20 mm and of different materials: (T) polyamide, (•) Teflon, 

(■) glass, (♦) steel and (a) tungsten carbide. ( ) Power 

law fit {5/d){H/d)-°-^ oc {p/pgf with /3 = 0.61. 



III. INFLUENCE OF LATERAL 
CONFINEMENT 



II. THE UNBOUNDED CASE 

We first present briefly our results obtained in the 
usual unconfined case: the container is large enough 
{D = 190 mm ^ d) so that the sphere is not affected by 
the surrounding walls and also high enough (6 = 300 mm 
^ S) to avoid any bottom wall effect as we shall see in 
the next sections. The influence of the different parame- 
ters on the penetration depth is shown in Fig. [5] for five 
different sphere materials with a density ranging from 
1.14 X 10^ (polyamide) to 14.97 x 10^ kg m~^ (tungsten 
carbide) and different sphere diameters d ranging from 
5 to 40 mm. Note that each point corresponds to an 
average of about ten experiments, and error bars cor- 
responding to the standard deviation are displayed only 
when larger than the symbol size. The penetration depth 
is larger for larger falling distance and larger projectile 
density. The data are well fitted by power laws of the 
form 



d \Pn J \ d 



(1) 



with A ^ 0.37±0.01, (3 = 0.61±0.02 and a = 0.40±0.04. 
The values of a and /3 are in close agreement with the 
previous studies already mentioned [a, ■ 



We are now interested in the influence of a lateral con- 
finement on the penetration depth, using cylindrical con- 
tainers of large height (6 — 300 mm) but with smaller 
diameters D ranging from 128 down to 24 mm. The pen- 
etration depth 6 is displayed in Fig. [3] as a function of 
the total falling distance H, for a steel projectile of di- 
ameter d = 19 mm falling in the different vessels. The 
two data sets corresponding to the two largest contain- 
ers {D — 128 mm and D — 190 mm) coincide, suggesting 
that for large enough vessel diameters the influence of 
the surrounding walls becomes negligible and the medium 
can be considered as unbounded in the radial direction. 
For smaller diameters D, data still fall into straight lines 
in the log-log plot of Fig. [3] but with smaller slope val- 
ues suggesting that Eq. ([1]) remains valid with smaller a 
values. The A and a values extracted from the fit of the 
data of Fig. [3] by Eq. (P) are reported in Figs. [3] and O 

The inset of Fig. 2] displays a as a function of D for 
the steel sphere of diameter d = 19 mm (•). A constant 
plateau value appears at large enough diameters D {D > 
100 mm) with the constant value a ~ 0.4. At smaller 
D, a decreases significantly with D and seems to tend 
towards zero when D approaches the sphere diameter 
d = 19 mm, as no deep penetration (larger than d/2) 
would be possible for D < d. Data for different sphere 
diameters do not collapse in the a vs D plot, but they do 
collapse in the a vs {D/d — 1) plot, as shown in Fig. H) 
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FIG. 3: Normalized penetration depth 5/d for a steel sphere 
(d = 19 mm) as a function of the normalized total falling dis- 
tance H/d for different cylindrical vessel diameters; (t) D = 
24 mm, (v) D = 35 mm, (■) D = 40 mm, (□) D = 50 mm, 
(♦) D = 62 mm, (o) D = 80 mm, (a) D = 128 mm, and 
(A) D = 190 mm. 
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FIG. 4: Power exponent a of Eq. ([T| as a function of [D/d—l) 
for steel sphere of diameter (v) d = 5 mm, (•) d = 19 mm and 
(□) d = 40 mm, and (o) tungsten carbide sphere of diameter 
d = 40 mm. (— ) Best fit by Eq. ^ with ~ 0.39 and 
\c, — 0.8. Inset: same data as a function of D. 



In this plot, data are well fitted by the exponential law 

D-d 



1 — exp — 



Xnd 
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where ctoo — 0.4 corresponds to the "infinite" (un- 
bounded) case and — 0.8 characterizes the range 
of wall effects. The characteristic distance Xad of lat- 
eral wall effects is thus found here a little smaller than 
one sphere diameter, which is rather small. Hence, for 
D/c? < 5, the surrounding walls play a key part and pre- 
vent a deep penetration of the projectile. For D/d > 5, 
the container has a vanishing influence and the sphere 
reaches a limiting penetration depth independent of D. 

The prefactor A follows the same kind of decrease for 
decreasing container diameters D (see Fig. [SJ, but we 



FIG. 5: Prefactor A of Eq. ((T| as a function of D/d with the 
same notations as in Fig. |4] ( — ) Best fit by Eq. (|3]) with 
Aoo — 0.37 and Xa — 0.7. Inset: same data as a function of 
D. 



believe that the relevant parameter for the lateral con- 
finement is here D/d rather than {D/d — 1), as A must 
vanish for vanishing D/d. The A variation is quite well 
fitted by the exponential law 



A = A., 



1 — exp 



D 

XAd 



(3) 



where A^o — 0.37 corresponds to the unbounded case and 
Xa — 0.7 characterizes the range of wall effects. Note 
that the Aq and A^ values are about the same as these 
two parameters characterize the same wall effect. 

Besides, the power exponent (3 for the density ratio 
[Eq. ([1])] does not depend on the lateral confinement but 
may be considered as constant, as data points for tung- 
sten carbide (o) twice denser than steel fall on the same 
curve in both Figs. 4 and 5. 

The reduction of the penetration depth by a lateral 
confinement is clearly due to an enhanced blocking effect 
by the walls (less radial dilatancy). The precise mech- 
anisms responsible for this reduction remain to be un- 
derstood. For instance, the non-linear "pressure" evolu- 
tion inside the granular packing due to the well-known 
Janssen effect, and valid both in static and dynamic sit- 
uations fails in explaining an enhanced force on the 
sphere in the confined case. 



IV. INFLUENCE OF NORMAL CONFINEMENT 

Let us now discuss the influence of a normal confine- 
ment by a close bottom wall. We investigate this ef- 
fect by varying the thickness b of the granular layer con- 
tained in a vessel of large enough diameter D to avoid 
the lateral wall effects discussed before. As we are in- 
terested in the sphere penetration inside the granular 
material, we restrict however our study to thick enough 
layers with b > d. The normalized penetration depth 
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FIG. 6: Normalized penetration depth 5/d of a steel sphere 
(d = 19 mm) as a function of the normalized total drop height 
Hjd, for different heights h of the granular layer: (v) hjd — 
1.2, (□) b/d = 1.9, (o) b/d = 2.6, (a) b/d = 15.8. Each 
horizontal dashed line corresponds to each initial height b of 
the grarmlar layer. ( — ) Power law fit {5/d) oc (H/d)" with 
a — 0.4, for the unbounded case {b/d = 15.8). 

S/d oi the steel sphere of diameter d = 19 mm is dis- 
played in Fig. [S] as a function of the normalized total 
drop distance H/d for three layer heights b ranging from 
23 to 50 mm (1.2 < b/d < 2.6) together with the mi- 
bounded case 6 — 300 mm {b/d = 15.8). The penetration 
depth remains unchanged by the presence of the bot- 
tom wall until the projectile approaches at a very short 
distance from the wall: 6 follows the unbounded curve 
until the sphere impacts the bottom wall indicated by 
the dashed lines in Fig. [B) This very short range effect 
is consistent with recent measurements of the force on 
a flat plate approaching in a quasi-static way the solid 
bottom boundary of a granular sample, indicating that 
the penetration resistance increases exponentially near 
the boundary [l^, At a high enough falling height, 
the sphere impacts the wall (a characteristic shock sound 
can be heard) and bounces back leaving a crater below. 
For thin enough layers as already studied in Ref. 
the opened crater does not get enough filled by inwards 



avalanches so that no grains are at the center and the fi- 
nal rest position of the sphere corresponds to 6 = b as the 
sphere comes to rest on the bottom "dry" wall [see (v) in 
Fig. [5]. For thicker granular layers, inwards avalanches 
are large enough to bring back grains down to the crater 
center so that the sphere comes back at rest at S < b 
[see (□) in Fig. [B]. In that case, the larger is the falling 
height, the longer is the rebound time and thus the time 
for avalanche, so that the thicker is the layer of back 
grains at the crater center: this explains the decreasing 
value of 6 for an increasing falling height in this part of 
the curve. Data for b/d — 2.6 (o) would also display such 
a behavior for larger falling heights. 



V. CONCLUSION 

We have measured the penetration depth 6 of spheres 
of diameter d impacting a granular medium contained 
in a cylindrical vessel of diameter D. This penetration 
depth may be strongly affected by the presence of sur- 
rounding walls. The presence of a bottom wall perpen- 
dicular to the sphere motion has a very short range infiu- 
ence on S as the sphere is only affected when the distance 
from the wall has reached a few grain diameters. On the 
other hand, the presence of lateral walls parallel to the 
sphere motion has a larger range influence with a typ- 
ical characteristic distance A ~ d. For D/d < 5, the 
walls have a strong influence and reduce significantly the 
penetration depth. For D/d > 5, they have a vanishing 
influence. The two (bottom and lateral) wall effects have 
been studied independently but it could be a hard cou- 
pling between them for geometrically strongly hindered 
grain motions. 
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